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Introduction: 

Combustion instability in solid rocket motors and liquid engines is a complication 
that continues to plague designers and engineers. Many rocket systems experience 
violent fluctuations in pressure, velocity, and temperature originating from the 
complex interactions between the combustion process and gas dynamics. During 
sever cases of combustion instability fluctuation amplitudes can reach values equal 
to or greater than the average chamber pressure. Large amplitude oscillations lead 
to damaged injectors, loss of rocket performance, damaged payloads, and in some 
cases breach of case/loss of mission. Historic difficulties in modeling and predicting 
combustion instability has reduced most rocket systems experiencing instability 
into a costly fix through testing paradigm or to scrap the system entirely. 
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Figure 1. Injector faceplate of the FI engine that powered the Saturn V rocket required over 1300 test firings to 
configure the injector baffles and reduce the pressure oscillations to 10% of the chamber pressure. 


During the early development of rocket propulsion technology scientist and 
engineers were cued to the underlying physics at play through the measurement of 
vibrating test stands, observation of fluctuating exhaust plumes, and most notable 
the audible tones that accompanied instabilities. These observations lead the 
pioneers of combustion instability research to acutely focus their modeling efforts 




on the acoustic waves inside combustion chambers. This focus on acoustics is quite 
logical given that the measured frequency of oscillation often closely matches the 
normal acoustic modes of the combustion chamber. But, this narrow focus misses 
contributions made by rotational and thermal waves that are a direct result of, or 
closely couple with, the acoustic wave. A more complete depiction of combustion 
instability oscillations is achieved when a global energy based assessment is used. 
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Figure 2. Pressure trace of a stable (red) and unstable (blue) solid rocket motor [1]. 


Recent advances in energy based modeling of combustion instabilities require an 
accurate determination of acoustic frequencies and mode shapes. Of particular 
interest is the acoustic mean flow interactions within the converging section of a 
rocket nozzle, where gradients of pressure, density, and velocity become large. The 
expulsion of unsteady energy through the nozzle of a rocket is identified as the 
predominate source of acoustic damping for most rocket systems. Recently, an 
approach to address nozzle damping with mean flow effects was implemented by 
French [2]. This new approach extends the work originated by Sigman and Zinn [3] 
by solving the acoustic velocity potential equation (AVPE) formulated by perturbing 
the Euler equations [4]. Determining eigenvalues of the AVPE, 
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is considerably more complex than the traditionally used pressure 
equation, 
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and requires numerical approximations of the chamber flow- field and eigenvalues. 


Modeling chamber gas dynamics: 

The latest theoretical models for oscillatory disturbances in high speed flows 
require a precise determination of the chamber acoustic eigenmodes. But first, a 
simulation of the mean flow properties of the combustion chamber must be 


performed. COMSOL multiphysics provides a numerical platform for conveniently 
and accurately simulating both the chamber gas dynamics and internal acoustics. 
This finite element software package provides multiple predefined physics modules 
along with a generalized mathematics module. 

The present study employs COMSOL multiphysics finite element framework to 
model the steady flow-field parameters of a generic liquid engine using the laminar 
flow equations within the High Mach Number Flow (HMNF) physics module. The 
HMNF module makes use of these fully compressible Navier-Stokes equation for an 
ideal gas, which are detailed in the COMSOL documentation. In order to account for 
the injection of hot gas due to the burning propellant, the injector face plate is 
modeled with a uniform inward flow of combusted propellant gas. All other solid 
boundaries are modeled with the slip boundary condition, and the exit plane is 
modeled with the hybrid outflow condition. 



Figure 3. Simulated liquid engine geometry with boundary conditions. 

Results from the HMNF analysis are reviewed to assure a valid and converged 
solution. Mean flow parameters such as pressure, density, velocity, and speed of 
sound are needed to model the AVPE. Of considerable interest are the values of the 
mean flow in the converging section of the nozzle, near the sonic choke plane. The 
sonic plane, where Mach number is equal to one, creates an acoustic barrier in the 
flow. In order to create an accurate geometry for the acoustic analysis the sonic 
plane (pictured in magenta) is extracted from the HMNF analysis. 



Figure 4. Velocity stream lines plotted over chamber pressure. The Mach 1 surface is plotted in magenta. 


Modeling chamber acoustics: 

The Coefficient Form Partial Differential Equation module of COMSOL multiphysics 
is used to determine the complex eigenvalues of the AVPE. Mean flow terms in the 
AVPE are supplied by the solution from the HMNF analysis. Gas dynamics within the 
combustion chamber play a key role in defining the boundary conditions for the 
acoustic analysis. Within the converging and diverging section of the rocket nozzle 
gradients of chamber pressure, velocity, and density grown becoming theoretically 
infinite at the sonic (Mach number equal to 1) plane. Downstream of the sonic plane 
acoustic disturbances are convected with the mean flow at speeds greater than the 
speed of sound. This condition prevents disturbances downstream of the sonic 
plane from propagating back upstream. The diverging section of the nozzle is 
acoustically silent, and does not affect the chamber acoustics. The simulation 
geometry is truncated at the nozzle sonic line, where a zero flux boundary condition 
is self-satisfying. The remaining boundaries are modeled with a zero flux boundary 
condition, assuming zero acoustic absorption on all surfaces. 
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Figure 5. Acoustic analysis geometry with boundary conditions. 


The eigenvalue solver produces complex eigenmodes and eigenvalues representing 
each acoustic mode and its complex conjugate. The real part of the complex 
eigenvalue represents the temporal damping of the acoustic mode with the 


imaginary part defining the frequency of oscillation. The complex eigenvectors 
represent the spatial amplitude and phasing of the acoustic wave. Comparing the 
acoustic mode shapes derived using the classic homogeneous wave equation 
(Helmholtz equation) to those derived using the AVPE demonstrates the benefits of 
higher fidelity models that correctly represent the underlying physics. Including 
mean flow terms in the AVPE accurately models the phase shift caused by the steady 
gas flow. Phasing is extremely important since combustion instability models make 
use of temporal and spatial integration of the acoustic eignvectors. 



Figure 6. Comparison of the first tangential Eigenmode calculated using the classic homogeneous wave equation 
(left) and the AVPE (right) of a half period (T) of oscillation. 


Utilizing COMSOL mutliphysics to simulate the rocket gas dynamics and acoustic 
eigenmodes provides a more accurate mode shape over the previously techniques. 
The higher fidelity acoustic representation is easily incorporated into combustion 
instability models to give rocket designers and engineers greater predictive 
capability. The inclusion of damping devices, such as baffle, or changes in operating 
conditions can now be more accurately modeled before testing. 

Continued work: 

A more complete depiction of combustion instability includes rotational oscillations 
and thermal oscillations in conjunction with chamber acoustics. Rotational 
oscillations occur as a direct result of the acoustic oscillation, where thermal waves 
can also be present in the absence of acoustic fluctuation. Continued work using 
COMSOL multiphysics will focus on solving the viscous rotational wave that 
accompanies all acoustic oscillations. 
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